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Abstract. Starting from an elementary model and refining it to take into account 
more realistic effects, we discuss the limitations and advantages of matter-wave 
interferometry in different configurations. We focus on the possibility to apply this 
approach to scenarios involving antimatter, such as positrons and positronium atoms. 
In particular, we investigate the Talbot-Lau interferometer with material gratings and 
discuss in details the results in view of the possible experimental verification. 


1. Introduction 

Matter-wave interference is at the heart of the quantum mechanical nature of particles. 
While this phenomenon has been observed for electrons [1, 2, 3], neutrons [4, 5], atoms 
and molecules [6, 7, 8] using a variety of different experimental tools, no experimental 
tests exist on elementary antimatter particles, or matter-antimatter systems. However, 
beams of antiparticles at low energy are becoming increasingly available, as in the case 
of antiprotons at the CERN Antiproton Decelerator [9] or in the case of positrons (and 
the associated positronium production) in 22 Na source systems coupled with Surko traps 
[ 10 ], 

In this paper we discuss the optical analogy and the main principles of Fraunhofer 
and Talbot matter-wave interference regimes, considering material gratings, in order to 
introduce the issues and the problems of antimatter interferometry. Positrons (e + ) are 
proposed as our Erst antimatter system to study and positronium (Ps) is the atom that 
we will be considering as a matter-antimatter symmetric system. The antiproton (p) 
case will also be shortly discussed. 

The paper is structured as follows. In section 2 we review the basic elements of 
quantum diffraction theory of particles from a grating and describe the build-up of 
the statistical interference pattern. Section 3 focuses on the incoherence due to the 
source, such as the effect of the particle velocity spectrum and the source geometrical 
extension. In section 4 we address the interaction between particles and a grating 
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Figure 1 . (Color online) A single particle of momentum p = mv impinging on an 
TV-slit grating (G); detection will take place on a screen, placed at a distance ( L ). The 
grating has a period D and the width of each slit is a. The 2 -axis is orthogonal to the 
x-y plane. 


considering both neutral particles and charged particles. Section 5 is devoted to Talbot- 
Lau interferometry: we describe the geometry, and its advantages with respect to single 
grating setups. Furthermore, we numerically show how the fringe visibility is affected by 
the particle velocity spread, when realistic parameters are used to carry out Monte Carlo 
simulated experiments. Finally, we close the paper drawing some concluding remarks 
in section 6. 


2. Basic quantum model of diffraction 


In this section we review the basics of matter-wave interferometry. We assume that a 
particle moving along the y- axis with de Broglie wavelength A = h/(mv), m and v being 
respectively its mass and its velocity along the y axis, interacts with an iV-slit grating 
laying in the x-z plane (see Fig. 1). Upon assuming that the slits are sufficiently large 
along the z-axis, so that diffraction is negligible along that direction, we can represent 
the state just after the grating at time t = 0 and y = 0, being t = y/v, as the following 
superposition state [11, 12, 13]: 

N 

lf( N \x,t = 0) OC ^2^ n (x,t = 0), (1) 

n= 1 


where if n (x,t = 0), n = 1,..., N, are the wave functions describing the particle passed 
through the n-th slit (we assume, as usual, that the slits are independent). For a system 
of identical slits with period D , we can write if n (x, 0) = ifo(x — nD). Indeed, the actual 
expression of if n (x,t = 0) is dictated by the characteristics of the diffraction grating 
and its interaction with the incoming particle. As the grating prepares the system in 
the state of Eq. (1), we can assume that the motion along the x-axis is governed by the 
free Hamiltonian: 


FT - Px 
H ’ S ~ 2m 


( 2 ) 


Therefore, the evolved state ii>(x,t ) is obtained by solving the Schrodinger equation 
with the Hamiltonian (2). In particular, the particle probability density distribution 
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along the x-axis on the screen at position y — L (the interference pattern), is given by 
I(x) = |= L/v )| 2 , where: 

1 


if( N \x,t) = 


f+00 

’ ^ ( /\2 

/ exp 

1 —00 



^ N \x/ 0) da/, 


(3) 


that is formally identical to the Fresnel integral of classical optics [14]. The most common 
approach found in literature [12, 13, 15, 16] is to adopt an “effective” point of view and 
postulate a convenient form for the initial single-slit wave function, for example: 

0) = a~ l X[-± + nD^+nD)(x) 
where a is the slit width and 


(4) 


Xn{x) = 


1 if a; G hi, 

0 otherwise. 


This is the quantum mechanical analog of assuming uniform illumination in the 
treatment of light diffraction. Another useful choice for the initial single-slit wave 
function is a Gaussian function centered on the slit interval with a suitable variance 
<7, namely (we drop the overall normalization constants); 

(x — nD ) 2 


ip n (x, 0) = exp 


4a 2 


(5) 


This choice is more convenient, as many calculations can be easily carried out 
analytically on Gaussian functions. In this case, the parameter er is usually set to 
a = a/{2\f2/n). Upon introducing the rescaled variables: 

ht LX 

4vra 2 ’ ^ ^ 

and considering a two-slit setup, the time evolved wave function outgoing from a double 
slit setup reads 


x = -, D = — , and L = 

cr a 2ma z 


fp2{x,L) = Cn.exp 
n= 1,2 


[X — X r 


4 1 + L 2 


1 - i L 


where X\ = —D/2 and X 2 = +D/2. The generalization for a set of N equally separated 
slits is straightforward. We have introduced the relative normalization constants C n , 
to account for a possible asymmetry in the beam preparation [13]. In the following we 
assume perfect symmetry C\ — C 2 — 1. After simple algebraic manipulations, defining 

(x ± D/2) 2 ' 


F± = exp 
the intensity reads 


2(1+ L 2 ) 


I(x, L) = F + + F_ + 2/FVFL 


cos 


LxD 


2(1 + L 2 ) 


(7) 
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Figure 2. The evolution of the interference pattern (7) shown with its dependence on 
the screen distance L. 


which clearly shows the appearance of an interference pattern due to the oscillating term 
(see Fig. 2). It is worth noting that the condition for observing the interference maxima 
in far field turns out to be the usual relation of classical optics, thus in the limit L 1, 
the condition for observing a maximum reduces to 


xD xDtt 

—— = 2nir —> - = 'U7T 

2 L LX 

which is indeed the expected classical relation. The formal analogy with classical optics 
[see Eq. (3)] also ensures that the choice of the initial single slit profile impacts only 
the envelope of the intensity pattern and not its oscillatory behavior. The classical 
Fraunhofer field outgoing a double slit setup reads 


/class (x,L) oc sine 2 



1 + cos 



while starting from (7) it easy to recover a Fraunhofer-like expression by taking the far 
held limit in the form L> 1 and x ^ D, so that 

/ £2 

F + = F_ ~ exp 


2 L 2 

and finally, in order to highlight the similarity with the classical expression in the 
Fraunhofer limit, we use Eqs. (6) in Eq. (7), obtaining: 


I(x, L ) = 2 exp 


- 2l2ira fi) 


1 + cos 






A Jj 


( 8 ) 


So, much alike the classical case, in a quantum treatment based on the free evolution 
of single-slit wave functions, the latter factorizes and determines the envelope of the 
pattern. 


3. Incoherence due to the source 

In order to describe a real experiment, the model introduced so far is not enough, since 
many relevant departures from the ideal situation arise. For instance, the particles’ 
speeds vary according to a given distribution, the particle source has a finite size and the 
collimation stage unavoidably introduces transverse momenta. Furthermore, focusing 
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Figure 3. Sketch of an interferometer operating in the far field, where an incoherent 
extended source (transverse size a s ) illuminates an fV-slit grating G (period D and slit 
width a) from a distance y s . 


on the scenario we are interested in, unstable antimatter atoms like Ps can decay in 
flight. All of these issues lead to incoherence effects. In general, if q = (qi, < 72 , • • •) is the 
vector of the physical parameters qi~ which can classically fluctuate in a real experiment, 
we can describe the overall incoherence effect by averaging the ideal intensity J(cc,f|q) 
given a suitable distribution p(q), that is: 


I(x,t) 


l /(x,t|q)p(q) dq 


(9) 


There are two relevant examples of incoherence: the one due to a finite transverse 
coherence length, the other due to the presence of a non-monochromatic beam. I 11 this 
section we focus on the first one, whereas the effects of a non-monochromatic beam will 
be considered in section 5 in the context of the Talbot-Lau interferometry. 

The experimental results show that the patterns of matter-wave experiments with 
multi-slit gratings can be described by considering a limited number of slits [15, 17]. 
We define an experimental parameter Z 0 , the coherence length , as the typical transverse 
length scale on the plane of the grating that sets how many slits can coherently take part 
to the interference process. From the physical point of view, a finite coherence length is a 
result of both the spatial extension and the intrinsic incoherence of the sources typically 
employed in matter-wave experiments. In order to take into account this effect, we add 
a (common) random transverse momentum k x along x-direction to the wave function 
ijjnfxjt = 0) associated with each slit, namely, ijj n (x,t = 0) exp (ix/c-c). If we assume 
that k x is distributed according to a Gaussian distribution with zero mean and variance 
1 //q, one obtains the following analytical result for the intensity, valid in the Fraunhofer 
limit [17]: 


l'{x,L) 

where 


2nN 

~JT 


if(x) 


N ~ V AT 

— n 

—^ 7 — exp 


n=l 


N 


(nDf 

~w 


cos 


/ 2nnDx\ 

v ) r 


f)(x) = / i/j(x') exp [ i-^- x' ) da/, 


( 10 ) 


( 11 ) 
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is an envelope corresponding to the rescaled Fourier transform of the single-slit wave 
function if(x) = ^Q{x,t = 0). Since for nD 3> Iq, the corresponding exponential term 
suppresses the interference, Iq can be regarded as the coherence length predicted in 
Ref. [17], which is inversely proportional to the transverse momentum spread. This 
could allow to determine the coherence length a priori, without resorting to a fit of 
the model to experimental data, ffowever the transverse momentum distribution is not 
always easily guessed. For example, if two successive slits are used as collimators, a 
bound on the maximum transverse momentum could be established with a geometrical 
construction [17]. Nevertheless, since we are interested in more compact geometries as 
in the presence of in-flight decay of unstable antimatter, limiting the dimensions of the 
apparatus will be of the utmost importance. So we would need an estimate on the 
coherence length for an apparatus of the kind of Fig. 3. 

We can obtain the intensity within the framework of the model given in Eq. (9) as 
follows. We assume that at random time a particle is emitted with a speed v from a point 
x s of the source (located at the distance y s from the grating), following a distribution 
p{x s ,v) that is determined by the nature of the source itself. After its emission the 
particle crosses the grating and produces an interference intensity pattern that depends 
parametrically on these quantities. Under the same assumption there discussed, the 
overall intensity at the screen is thus given by Eq. (9), that now reads: 

I(x,L) = ! I(x, L\x s ,v) p(x s ,v) dx s dv. (12) 


The integration can be performed via Monte Carlo (MC) method, as it scales well 
with the dimension of the parameter space. Moreover, we can refine our analysis, 
e.g., taking into account the instability of the particles and their lifetime. As a first 
approximation, we could simply discard the particles that do not reach the detector 
plane. This corresponds to employ a detector able to discriminate between a true event 
and the background noise induced by the decay in flight. 

In order to obtain the same results as in Eq. (10) for a suitable choice of l 0 
and in the Fraunhofer approximation, we should consider a monochromatic beam and 
average only over the source dimension x s , assuming a uniform probability density 
p(x s ) = af 1 X[-a s / 2 ,a s / 2 ]{x), a s being the source dimension. It is worth noting that in 
our simulations the average intensities are computed retaining the full accuracy of the 
Fresnel integral, i.e., without the Fraunhofer approximation. In a setup like the one 
shown in Fig. 3 a comparison between Eqs. (10) and (12) shows that the coherence 
length l 0 can be estimated as [18]: 


(13) 
2(7 s 

The dependence on the physical parameters is in agreement with the naive estimate 
associating the coherent illumination region in this kind of setup with the width of the 
central diffraction peak for a slit of size a s , where a s is the transverse extension of the 
source [18]. In Fig. 4 we show the Monte Carlo simulations of the far held interference 
pattern for different values of the source dimension a s and a particular choice of the 
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Figure 4. Monte Carlo simulation of the far field setup (Fig. 3) from Eq. (12) for a 
monochromatic beam of Ps atoms with A = 3.6nm (mp s = 2 m e , v = 10 5 m/s) and 
different values of the source dimension er s = 900,90 and 9 /xm, corresponding to the 
coherence length Iq = 1,10 and 100/xm (as shown in the plot). We also set D = 10/im, 
a = 3 ym, N = 10, y s = 0.5 m and L = lm. Interference disappears when Iq < D , 
and the contrast starts to decrease when Iq ss D. The dashed curves refers to the 
corresponding single-slit diffraction envelopes. 


other involved parameters. For the simulations we considered a typical Ps velocity 
v = 10 5 m/s [19], which leads to A = 3.6nm. We note that as the coherence length 
approaches the critical value D, we observe a decrease in the contrast or visibility of the 
pattern (see Fig. 4): 


C 


hi lax f min 
hri ax 4“ A n i n 


(14) 


defined as the difference between the intensity of a maximum and its adjacent minimum. 
We will reconsider the implication for the design of an experiment in sect. 5. 


4. Interaction with material gratings 

In the previous section we addressed the interferometry problem assuming that the 
particle did not interact with the grating. As physical quantum mechanical objects the 
particles interact in various ways with the walls of the material grating. The formalism 
we developed so far is sufficiently general to account for this effect by modifying the 
initial wave function accordingly. 

In the previous sections we have seen that the fundamental building block for 
quantum models of diffraction from a grating is the single slit outgoing wave function 
0), which is usually postulated to be of either Gaussian or rectangular shape. If 
the potential V (£, y) acting on the particle is known in the region within one slit, we can 
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Figure 5. Trapezoidal bars cross section with the narrower side of the slit facing the 
beam, displayed for a wedge angle /3 and nominal maximum and minimum widths om 
and a m = clm — 26 tan /?, respectively. 


account for the interaction by treating the slit, as a phase mask , producing a transmission 
function of the form £a(0 = e 1<p ^) [20, 21, 22], implying the substitution (see Fig. 5) 

Mt, 0)^^n(e,0)e i ^ ) . (15) 

The standard approach [20, 22, 23] is now to determine the phase shift <p(£) via the 
semiclassical eikonal approximation. Denoting with v the particle speed we can write: 

¥>(0 = -^ J V(£,y)dy. (16) 

As discussed, the Fourier transform of the single slit wave function sets the envelope 
of the diffraction pattern [see Eq. (11)], which in the absence of interactions depends 
on the nominal (geometrical) width of the slit a. We will show that as a first order 
approximation the effect of a potential is a reduction of the effective slit width; this has 
been observed in various situations [20], a notable example being the Cgo experiments 
[6]. We find: 

VKO ex P ( 17 ) 

where x denotes the transversal coordinate on the screen plane, consistently with our 
notation. As suggested by Ref. [20], we approximate the above Fourier transform using 
a cumulants expansion [24] [a normalization of if(£) ensuring ./’ "0(0 d£ = 1 is implied] 

log ip(x) = K n -—— with c = -ry- (18) 

n\ XL 

n= 1 

where the cumulants n n are defined in terms of the raw moments /q, of namely, 

W = J 

and we are interested in the first two terms only 
K\ = p>i, and K -2 = p -2 — i 


( 19 ) 
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which, truncating the series to the second order, result in the approximation 


( C 2 

ip(x) S3 exp I IKjC - k 2 — 


( 20 ) 


\ip(x)\ 2 S 3 exp {—2()3m[Ki] — Se[fv 2 ]C 2 } 


( 21 ) 


In order to give a physical meaning to the expansion parameters, we analyze the case 
of no interaction, with a rectangular wave function 0(£) = a$ 1 X[-a 0 / 2 ,a 0 / 2 ]- The first 
raw moment (the mean) of this distribution vanishes due to parity " 0(0 = 0 (—£)> an d 
in general we expect this to be true in any realistic situation, as both a reasonable 
interaction potential and the wall geometry will be symmetric with respect to £ —» — 
The second moment is simply evaluated 


k 2 — h 2 


1 

Qo 



o| 

12 


leading to: 




exp 


1 fTraox^ 2 

3 V XL ) 


As one may expect, the exact expression is the well known = sinc[ 7 raa;/(AL)], the 
central peak of the sine function coincides with the Gaussian approximation, suggesting 
to identify an effective slit width as follows: 


a efF = yi2!Re[K 2 ]. 


( 22 ) 


We will now distinguish the two cases of neutral and charged particles in term of the 
potential. 


4-1. Neutral particles 

The nonretarded van der Waals atom-surface potential [20], which affects all types of 
neutral polarizable particles, is expressed in terms of the distance from the surface and 
the coefficient C 3 as 

Kdw(r) = -^, (r > 10 A). ( 23 ) 

It has been shown with direct electron microscope imaging, that diffraction gratings 
of the type commonly used for matter-wave experiments can have a trapezoidal (see 
Fig. 5) slit profile [20, 22, 23, 25, 26], as a result of the fabrication process. Therefore 
it is useful to study this immediate generalization, from which the trivial parallel-plane 
profile is recovered in the limit {3 —> 0. Calculations also show that the introduction 
of even a small wedge angle has a significant impact on the effective width of the 
slits compared to parallel-planes approximation at the same thickness and interaction 
strength. 

The potential can be written in terms of the distance of the generic point (£, y ) 
from the right and left grating walls, respectively d\ and d 2 . The symmetry of 
the system implies that these quantities are related by d 2 {f 1 y) = d\{— £, y). Since 
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di(f, y) = a m /2 — £ + y tan /3, as it is evident from Fig. 5, the projection of the distance 
on the normal vector to the side wall is simply obtained by multiplying d 2 (£,,y) for the 
cosine of the wedge angle (3. The integration in Eq. (16) is straightforward, and yields 
for the phase shift of the potential (23) 


7h (£,/?) 


Cs 

v 


2 {h cos 3 (3 tan (3) 
a m - 2 (£ - ytan/d)] 2 


y=S 


y =o 


Taking both surfaces into account we obtain: 


(24) 


= (25) 

We note that at fixed geometry the ratio R = C '3 / v sets the overall scale of the 
interaction strength, as tp(f,fd) oc C 3 /v. Numerical estimates of Eq. (22) with the 
phase shift (25) are in good agreement with the experimental and theoretical results in 
[20]. Therefore, this kind of interactions can be accounted for by a reduction in effective 
slit width, at least up to the highest interaction strength tested experimentally, that is 
R max = C 3 /v ~ 2.74 • 10 -12 meV • nm 2 (this is obtained for Kr atoms at v = 400 m/s 
with a SiN-j, grating, and a measured C 3 = 1.1 meV ■ nm 3 ). Hoinkes’ empirical rule [27], 
which has been confirmed experimentally [20], states that for a given material C 3 is linear 
in the particle static polarizability a. The static polarizability of Ps atoms, a/ Ps \ in 
the ground state (estimated treating it as an hydrogen-like atom with the appropriate 
reduced mass) is cA Ps ) = 8 cA H ) ~ 5.33 A , roughly twice that of the Kr atoms. In 
turn, Ps atoms have an interaction scale R < R max down to speeds of ~ 800 m/s, 
corresponding to a very low energy E ~ 3.6 -10 -3 meV, while the lower speed limit for an 
experiment with antihydrogen atoms (assuming = a^) would be lOOm/s. 

Therefore, we can safely conclude that a treatment of the van der Waals interaction in 
terms of Eq. (22) is fully adequate to describe experiments involving current Ps sources 
[28, 29], 


4-2. Charged particles 

The above procedure is easily generalized to all potentials depending on the distance 
from the surface as oc r~ n . Relevant examples are the retarded van der Waals interaction 
(n = 4) and the electrostatic potential (n = 1). In the AEgIS experiment at CERN 
[30], the production of antihydrogen and Ps atoms also involves, as an intermediate 
step, the realization of a steady beam of charged antimatter, specifically e + and p. 
It will be interesting to carry out interferometry experiments on these objects as well, 
because no successful demonstration of interference has been obtained for these systems, 
yet. Moreover, to the best of our knowledge, this statement also applies to any kind of 
charged system heavier than an electron. 

It is a standard result in electrostatics that the potential acting on a point charge 
q sitting at an orthogonal distance r from a dielectric surface (relative permittivity e) 
is that of a point charge q' = g (1 — e) (1 + e) 1 located on the axis of symmetry with 
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respect to the plane of the surface. Therefore, for the same geometry represented in 
Fig. 5, we can argue that the potential will be given by 


VJr) = 


1 — e q 2 


1 + e 47re 0 2 r 

where r the distance from the grating wall. By using Eq. (16) we End: 


<p? (£,?) = 


q\ l~e)(l + e) 

87reo hv sin /3 


-l 


log 


2 (£ — <5 tan /3) 

Orn 2q 


(26) 


(27) 


where a m , 5 and f3 are the same as in Fig. 5. Once the total phase shift = 

/3) + /3) is obtained, as in Eq. (25), we can follow the same procedure and 

account for the electrostatic interaction introducing an effective slit width, given by 
(22). However, the orders of magnitude involved might be very different. If we compare 
the potential strengths at the center of a slit of width a = a m = (Im (for the sake of 
simplicity we set /3 — 0 ), then we obtain the following expression for the ratio: 

I4i q 2 1 — e a 2 
Kdw 47reol + eC3’ 

which shows a quadratic dependence on the slit width, descending from the different 
power-law scaling of the potentials. Assuming that q = e, the electron charge, we also 
have e 2 (47reo ) _1 = 1439.964 meV■ nm. Consistently, we recall from the previous section, 
that C 3 is of the order of a few meV ■ nm 3 , in turn we have V e \/ K-dW ~ a 2 10 3 [nm]. 
Though the difference seems very large, this is only a pointwise estimate, while the 
effective slit width is determined by the behavior of the phase shift over the whole range 
of £. Nevertheless, we can see that there are realistic situations where the calculated 
impact of electrostatic interactions is indeed very high, as discussed in table 1. We see 
that in view of the higher typical potential strength in comparison with the van der 
Waals interaction, for low energy antiprotons Eq. (22) predicts sizeable reduction in 
effective slit width, i.e., > 50%. It turns out that probably the interaction is too strong 
in this regime to be treated with this approximation. 


Energy [keV] 

a e fj e + nm] 

a e ff p nm] 

0.1 

401.3 

148.1 

1 

477.2 

285.8 

10 

497.1 

397.4 

100 

499.7 

460.0 


Table 1. Calculated effective width for realistic parameters (order of magnitude) 
applicable to possible experiments, ao = 0.5 /im, /3 = 5°, e = 4, for e + and p of varying 
energy. The grating thickness is set to S = 500 nm and <5 = 160 nm respectively, 
reflecting the typical scale necessary to absorb the particles completely outside the 
slits for silicon at 1 keV reference energy. 


Exact numerical calculation of the envelope function [still using the eikonal 
approximation for the phase shift (16)] shows that for low energy particles there is indeed 
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a stark departure from the sine shape expected for weak interactions. Experimental data 
with electrons in the 0.5 4-4keV energy range exist [22]: they refer to the same grating 
geometry described here, where although no explicit information on the effective slit 
width is reported, a model for an envelope function is developed and an analysis of the 
plotted data seems compatible with our theoretical predictions based on (27), which 
yield a small |o e ff — ao|/ao < 10% for that energy range. 


5. Talbot-Lau interferometry 


In section 3 we have discussed how the contrast of interference patterns is affected by the 
ratio between the coherence length / 0 and the grating period D. This imposes technical 
constraints on the design of an interferometer using the geometry shown in Fig. 3. First 
of all the finite resolution of the detector has to be taken into account. This parameter 
greatly depends on the kind of detector and on the particles involved. For example, in 
the case of anti-hydrogen, e + and p emulsion detectors could be employed, which are 
capable of a spatial accuracy up to 0.6 4-1 /jrn [30, 31]. Being L the grating to detector 
distance and D the grating period, the Fraunhofer diffraction orders are separated by 


Ax = 


LA 

~D' 


If Sx is the experimental sensitivity, in order to resolve each maximum of the diffraction 
pattern within at least an interval M 5x, with M integer, we should have 

^ >M8x^L>MD^, 

U A 


which imposes a constraint on L. It is clear that for a fixed wavelength and geometry 
both increasing M and reducing the distance L, which is of utmost importance with 
decaying particles, requires a decrease in the grating period D. 

Moreover, starting from Eq. (13) and requiring that the coherence length l 0 is at 
least M times the grating spacing D, we obtain the following condition on the source¬ 
grating distance: 


y s > MD 


2 <Js 

A ' 


Therefore, to obtain a good coherence either the distance y s has to be increased or the 
source dimension a s reduced as much as possible. Note that, apparently, by reducing the 
period D we can satisfy both conditions on L and y s , however a reasonable small value 
for D is fixed by the grating construction constraints. This poses technical challenges 
due to the particle decay in the first case (for ortho-Ps atoms with a lifetime r = 142 ns 
[32] and realistic thermal speeds v ~ 10 5 m/s [19], y s should be in the range of a few 
centimeters), and to difficulties in manipulating the beam size in the latter. In view 
of these consideration, we suggest that a different kind of interferometer would be best 
suited for experiments with antimatter, namely a Talbot-Lau setup [33, 34], which is 
sketched in Fig. 6. 
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Figure 6. Sketch of a Talbot-Lau interferometer. The first grating (Gi) is illuminated 
by an incoherent beam of mean wavelength A, and acts as an intensity mask providing 
the necessary coherence for illuminating the second grating (G 2 ). Their separation L 
is set to the observation distance from G 2 and the gratings have the same period D 
and slit width a. If L matches the Talbot Length If = D 2 / A, this setup produces on 
the detector plane high contrast fringes with period D. 


It is worth noting that the usual Talbot-Lau configuration involves a third grating 
as a scanning mask [33, 35], that is not necessary in our case since we assume that 
the high resolution of the detector will allow to directly resolve the diffraction pattern. 
The fundamental property of this geometry is that it produces high contrast fringes 
regardless of the coherence and spatial extension of the illuminating beam. There are 
two physical phenomena governing this apparatus: the Talbot self-imaging effect [34], 
stating that in the Fresnel region of a coherently illuminated periodic grating self-images 
of the grating transmission function will appear at L — uTl = nD 2 /X (as well as rescaled 
sub-images with a fractional period for half-integer multiples ), and the so called Lau 
effect [36]. This effect can be understood as arising from an incoherent superposition of 
patterns produced by laterally displaced, mutually independent, point sources, the role 
of which in the apparatus of Fig. 6 is played by the slits of the first grating [34], The 
periodic images thus produced can overlap “constructively” if the first grating has a 
suitable periodicity; under these conditions the elementary displacement on the screen 
plane produced by moving between adjacent sources equals the Talbot image period 
or an arbitrary integer multiple of the latter. In particular, this “resonance” condition 
is met in the configuration of Fig. 6 when L = T^. Geometrically, this setup bears a 
strong similarity to a classical moire deflectometer [37]. What discriminates between 
the purely classical and the quantum interference regime is the condition for diffraction 
to be negligible, namely [37]: 


LX 

~D~ 


< D. 


(28) 


A moire deflectometer and a Talbot-Lau interferometer as defined in Fig. 6 have 
in common that they produce a fringe pattern with period D. The question that now 
naturally arises is: how can the experimental results prove that the observed fringes are a 
true interference effect and not simple classical geometrical shadow patterns produced by 
ballistic particles? As mentioned, high contrast fringes are expected only if the grating 
separation is an integer multiple of the Talbot length, while for “classical projectiles” 
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Figure 7. Monte Carlo simulation of the fringe visibility modulation (contrast) as a 
function of L/Tl (Log scale), for a well defined particle velocity. For definiteness the 
following parameters (realistic for an experiment with 1 keV antiprotons) were chosen 
D\ = D 2 = D = 265 nm, N = 40 and a = 90 nm. The calculated Talbot length for this 
period and energy is T f = 77.9 mm. Note that the typical contrast peaks disappear 
when L/Tl <C 1: this corresponds to a classical behavior of the particles. 


the contrast does not depend on this condition. This property ultimately descends from 
the longitudinal periodicity of the so-called Talbot carpet , which is a distinctive feature 
of diffraction in the Fresnel region. Therefore, the observation of this kind of pattern is 
a proof of the wave character of the interfering particles. From the experimental point of 
view, this can be done by continuously adjusting the grating separation or changing the 
particle energy (hence the Talbot Length) in a monochromatic beam, and measuring the 
modulation in contrast as a function of the parameter L/Ti. If the apparatus is truly 
operating as an interferometer and not as a classical device, distinct peaks in contrast 
should be detected, as shown in Fig. 7 [38]. Recalling (28), we see that the classical 
limit corresponds to L C 7^, as it is confirmed by numerical calculations showing a 
weak dependence of the contrast on L in this region. This is clear from Fig. 7: as 
the grating distance L falls below Tl = 77.9 mm (given by the simulated period and 
energy) the contrast peaks disappear thus revealing a classical behavior of the particles. 
We note that in the simulation of Fig. 7 we assumed an incoherent particle beam, thus 
we applied Eq. (12) treating the slits of the first grating as a collection of extended 
incoherent sources. 

Another advantage of Talbot-Lau configuration is its robustness with respect to 
external stray fields. Suppose for example that the interferometer is subjected to uniform 
external electric and magnetic fields which exerts a force F on the charged particles: 
this force will be negligible if the corresponding deviation from a straight trajectory is 
smaller than the typical size of the finest structure in the observable pattern. Let us 
call this quantity A, and introduce the flight time of the particle r = L/v. Using the 
identity A = F crit r 2 /m, we obtain the following relation: 

_ hW 

-*crit 


mL 2 A 2 ’ 
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where for a Talbot-Lau geometry A — D. If we set L — Tl — D 2 / A, we have: 

Urn = -M, (29) 

mV 6 

that is independent of the particle energy. Furthermore, from the Lorentz force 
F = q(E + vB) one can deduce the critical values of the involved electric and magnetic 
fields by the simple relations E crit = F crit /q and B crit = F crit /qv. It is worth noting that 
in the case of Fraunhofer interferometry (far held) a similar calculation with A = LX/D 
leads to a critical force F^j = h 2 /(mLDX) oc v which is thus affected by the particle 
energy. 

Overall, the Talbot-Lau geometry has several advantages which are especially 
relevant for anti-matter interferometry. In particular, it allows to minimize the total 
length of the apparatus, a crucial requirement for decaying Ps, and to employ a larger 
source with weak coherence requirements, significantly increasing the particle flux. 

5.1. Realistic application 

We now present a realistic estimate of the expected contrast signal in an e + experiment, 
following the experimental methodology outlined in this section (see Figs. 6 and 
7). We assume that the particle mean energy Eq can be tuned between 5keV and 
20keV (reasonable for current continuous e + beams) with a narrow Gaussian energy 
distribution ( cr/E 0 < 2%). By varying the energy we scan the ratio L/T Ll where the 
Talbot length reads: 

D 2 y/2mEo . . 

T l = --. (30) 

A contrast peak is expected around L ps Tl- The above formula makes it clear that to 
scan this region there are specific complications related to each choice of which parameter 
to vary: the energy is bounded by technical constraints and furthermore provides a sub- 
linear scaling, whereas the grating distance can be varied arbitrarily. There are however 
technical complications in physically moving an apparatus sensitive to alignment over 
considerable lengths. 

We recall that in this configuration D sets the periodicity of the interference pattern, 
which should be larger than the detector resolution. However as evident from (30), the 
Talbot length scales quadratically with D. Thus it rapidly becomes very large, and too 
long an apparatus poses additional challenges related both to grating alignment and 
shielding of a larger region from stray fields. We set D = 2 /mi corresponding to a Talbot 
Length T L = 0.326 m at the median energy of the considered energy range, namely 
E 0 = lOkeV. Using Eq. (29) we obtain E CTit pc 0.2 V/m and H crit ~ 0.3 rnG for the 
maximum tolerable electric and magnetic held, respectively, evaluated at a Eq = 5 keV, 
which corresponds to the worst scenario. 

The maximum magnetic held is particularly critical, as the requirement is smaller 
than the natural magnetic held of the Earth, however considering that experiments with 
electrons and similar length scales involved have been successfully carried out [39], we 
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Particle energy [keV] 

Figure 8. Monte Carlo simulations of the expected contrast modulation as a function 
of E 0 , for the monochromatic case and with a Gaussian distribution of increasing a 
centered on E 0 . The total number of slits is set to N = 40, sufficiently large to provide 
realistic results. See text for more details on the parameters. 


believe that an appropriate mu-metal based shielding will be enough to circumvent the 
problem. Moreover an uniform constant magnetic field will only rigidly shift the pattern 
in space [40], in fact the above limits indicate the maximum allowed fluctuations (either 
in time or in space) of the E and B. 

Another problem our theoretical analysis allows to account for is the electrostatic 
interaction with the grating walls. First of all we have to set the slit width, and thus the 
open fraction of the grating; this is best set at a/D 30%, implying a = 0.6 pm. Higher 
values could improve the total particle flux minimizing the losses inside the material 
grating, but will also reduce the contrast due the overlapping diffraction peaks. We 
remark that for the gratings to work as true intensity masks, their thickness must be 
sufficient to stop all the positrons outside the slits. For such low energies few microns of 
SiN x will be sufficient. Therefore, assuming these parameters and a small wedge angle 
f3 = 10°, a grating thickness 5 = 800 nm and E 0 = 5keV, equation (22) predicts an 
effective slit width a e ff = 0.598 pm. This deviation is very small (« 0.5%), thus the 
effect will be completely negligible for the considered choice of parameters. 

Having evaluated the most relevant effects and set the geometry of the apparatus, 
we can now use Eq. (12) to predict the measured contrast modulation [38]. The result 
is shown in Fig. 8, which contains interesting indications. On the one hand we can 
see that a particle energy between 5 keV and 20 keV is enough to observe a full peak. 
On the other hand, we see that also a broader energy distribution (a = 0.25 keV or 
cr = 0.5 keV in Fig. 8) still allows for a good visibility of the contrast modulation. 
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6. Conclusion and outlook 

In this paper we reviewed the basic elements of diffraction theory applied to matter- 
wave interferometry. In particular, we focus our analysis on the possible issues arising 
from the use of charged or neutral antimatter particles. To sustain our investigation 
we also performed Monte Carlo simulated experiments based on realistic parameters. 
In particular, we have considered the effect due to a realistic source (extended and 
non-monochromatic) and the interaction with the grating as well as the influence of 
stray electromagnetic fields. We also found that van der Waals interactions with the 
material grating become critical for highly polarizable particle systems. In this scenario 
a possible solution could be resorting to light gratings [41, 42], We have shown that the 
better configuration to carry out matter-wave interferometry with decaying particles is 
given by the Talbot-Lau setup also in the presence of a Gaussian distribution of the 
particle energy, which realistically describes the actual e + and p beams. Furthermore, 
exploiting the high resolution capabilities of the antiparticle detectors, such as the 
nuclear emulsions, we have shown that the typical Talbot-Lau setup involving three 
gratings can be reduced to a two-grating configuration which indeed simplifies the 
experimental implementation. Our analysis paves the way to further investigations in 
order to design an experiment to demonstrate antimatter-wave interference also in view 
of possible applications in the emerging field of gravity experiments using antimatter 
[41, 43, 44], 
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